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Abstract Pulsed-field-gradient nuclear magnetic resonance (PFG-NMR) is used to obtain 
the true hydrodynamic size of complexes of peptides with sodium dodecyl sulfate SDS 
micelles. The peptide used in this study is a 19-residue antimicrobial peptide, GAD-2. Two 
smaller dipeptides, alanine-glycine (Ala-Gly) and tyrosine-leucine (Tyr-Leu), are used for 
comparison. We use PFG-NMR to simultaneously measure diffusion coefficients of both 
peptide and surfactant. These two inputs, as a function of SDS concentration, are then fit to 
a simple two species model that neglects hydrodynamic interactions between complexes. 
From this we obtain the fraction of free SDS, and the hydrodynamic size of complexes in a 
GAD-2-SDS system as a function of SDS concentration. These results are compared to those 
for smaller dipeptides and for peptide-free solutions. At low SDS concentrations ([SDS] < 
25 mM), the results self -consistently point to a GAD-2-SDS complex of fixed hydrodynamic 
size R =(5.5 ± 0.3) nm. At intermediate SDS concentrations (25 mM < [SDS] < 60 mM), the 
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18 apparent size of a GAD-2-SDS complex shows almost a factor of two increase without a 

19 significant change in surfactant-to-peptide ratio within a complex, most likely implying 

20 an increase in the number of peptides in a complex. For peptide-free solutions, the self- 

21 diffusion coefficients of SDS with and without buffer are significantly different at low SDS 

22 concentrations but merge above [SDS] =60 mM. This concentration is identified as an onset 

23 of crowding beyond which it is impossible, even in principle, to extract information about 

24 hydrodynamic size of the peptide-surfactant complex. 

25 Keywords Antimicrobial peptide • Peptide-micelle complexes ■ NMR diffusometry ■ NMR 

26 relaxometry 



27 Introduction 



Membrane-associated proteins and peptides are often studied in a micellar environment 
( Tulumello and Deber[ 2009 Sanders and Sonnichsen 2006| . Like membrane bilayers, mi- 
celles provide a hydrophobic-hydrophilic interface, but imlike them, they are small enough 
to enable solution NMR signals to be observed. Micelles are commonly employed in NMR 
iptructure determination of membrane proteins ( [Qureshi and Goto 2012 Tulumello and 
|Deber[[2009| , but have also been used in studies where the protern-lipid interaction itself is 



the focus ( [Cozzolino etall|2008|porein et al|[T996l|Yu et aIl[2006l|Romani et al||20T0) . NMR- 
based techniques have been utilized to study an important class of membrane-associated 
proteins that are called antimicrobial peptides (AMPs). 

AMPs are often short peptides consisting of 12 to 50 residues and act by interacting 
with (and often disrupting) membranes. AMPs have been shown to play an important 



role in attacking and killing microbes such as bacteria, viruses, and fimgi (Zasloff 2002 



Nicolas[|2009{|Hoskin and Ramamoorthy}|2008{|Chinchar et aH|2004| . Moreover, some AMPs 



exhibit activity against tumor cells in a mammal's body by disrupting the membrane of the 
diseased cells and targeting the cell interior without affecting the membrane of host cells 
( Rege et"aH 2007 1. This selectivity, for microbial and/or tumor cells, is thought to arise due 
to the amphiphilic structure of the AMP that has an affinity to the lipid bilayer structure of 
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45 the microbial cells as well as due to the interaction between the positive charge on the AMP 



with the anionic components of the tumor or pathogen cell membrane ( Epand and Vogel 



T999]|. Therefore, anionic sodium dodecyl sulfate SDS surfactant micelles are commonly 



employed in the structural studies of AMPs, as well as other membrane proteins < |Wang 



2008} [1999) [Whitehead etal] [2001) prfi et al} [1998) [Begotka etal] [2006) peaton et al} [2001 
Whitehead eraI)[2004)[Gao and Wong) [1998) [Buchko et al[[T998) . 



A knowledge of the hydrodynamic size of proteins plays an important role in un- 
derstanding their conformations ( [Jones et"al 1997| . This is also the case for peptides in 
peptide-micelle complexes, where there could be many coexisting conformations. The hy- 
drod5mamic size of complexes can be obtained by measuring diffusion coefficients and 
using the Stokes-Einstein-Sutherland equation Rh = KBT/67irjDo. This approach, however, 
is only strictly valid when the self-diffusion coefficient Dq is obtained by measuring the 
diffusion coefficient as a function of the surfactant concentration and then extrapolating to 
infinite dilution. Such a procedure is often not practical when the amount of peptide or pro- 
tein is limited in quantity. As a result of this, "apparent" hydrodynamic radii are routinely 
reported, without such extrapolation, in systems with rather large surfactant concentrations 
( [Binks et al)[T989}|Gim"el and Brown) [19%) [Sarker et al)[20Tl] . 

An important phenomenon to consider with respect to large macromolecular concen- 
trations is crowding. Macromolecular crowding refers to the non-specific excluded volume 
(steric) effect of macromolecules with respect to one another in an environment where the 



^acromolecular volume fraction is large; an example is a living cell with O=A0% (Zhou 
[et al[[2008| . At finite dilutions there are hydrodynamic corrections to diffusion (jBatchelor 



19761 even for a simple colloidal system of spherical particles. It is now being realized 



that electrostatic and hydrodynamic interactions sensitively affect macromolecular d5mam- 
ics ( [Zhou et al[|2008)[Schreiber et al[[2009) . 

The nature of the association of peptides with anionic SDS micelles depends on the 
details of the electrostatic environment; for example, cationic peptides bind more strongly 
than their zwitterionic coimterparts (Begotka et al 2006| . NMR diffusometry studies have 
found that peptide binding with anionic SDS micelles and zwitterionic dodecylphospho- 



4 



Suliman Barhoum, Valerie Booth and Anand Yethiraj 



74 choline (DPC) micelles are different, also due to the difference in electrostatic environment 

75 (Whitehead et al 2004[ |. Similarly, it was found that a cell-penetrating peptide (CPP) al- 



76 ters the dynamics and size of neutral and negatively charged bicelles in different ways 

77 ( [Andersson et al||2004) . 

78 PFG-NMR studies have shown that the hydrophobic interaction can play a signifi- 



ipant role on the binding of peptides and tripeptides to micelles (Deaton et al 2001 Orfi 



80 |et aH|1998} , as well as neuropeptides to a membrane-mimic environment ( [Chatterjee et al 



2004| . NMR studies were also carried out to explore the binding of a neuropeptide to 



82 SDS micelles in the presence of zwitterionic 3-[(3-cholamidopropyl) dimethylammonio]- 

83 1-propanesulfonate (CHAPS) surfactant as a crude model for cholesterol in the biological 

84 membrane. These studies showed that having comicelles composed of SDS and CHAPS sur- 

85 factants inhibits the hydrophobic interaction of the neuropeptide with the core of comicelles 

86 ( [Whitehead eFaIl[200T) . 



87 Since AMPs are subjects of much interest and also represent an even larger class of 

88 amphipathic, helical peptides, a helical AMP, GAD-2, was selected for this study. GAD-2 

89 and a related peptide, GAD-1, have been identified in recent efforts to discovery new AMPs 

90 ( [Fernandes et al 2010 [Browne eTal 2011 [Ruangsri et al 2012| . In this work, we used NMR 



91 diffusometry and relaxometry to study the interaction between the cationic GAD-2 AMP 

92 and an anionic SDS micelle as a membrane mimic environment. In order to do so, we use 

93 a simple mathematical model that is utilized to signal the changes in the nature of the 

94 macromolecular complexes in a system of nonionic polymer-anionic surfactant system in 

95 aqueous solution (Barhoum and Yethiraj 2010) . We compare the nature of the resulting 



96 peptide-SDS complex with those that form with two much smaller peptides, and are able 

97 to identify important distinguishing characteristics. We find, reassuringly, that the most 

98 minimal model to extract hydrod5mamic size works well for peptides, at least for those 

99 with the size (19 residues) of GAD-2; however, one must be careful to avoid the onset of 

100 crowding in order to reliably use these simple models. 
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101 1 Materials and Methods 

102 GAD-2 peptide with average molecular mass Mw=2168 g/mole was synthesized using solid 

103 phase chemical synthesis employing O-fluorenylmethoxycar- bonyl (Fmoc) chemistry, on 

104 a CS336X peptide S5mthesizer (C S Bio Company, Menlo Park, CA) following the man- 

105 ufacturer's instructions. The peptides were S5mthesized at a 0.2 mmol scale with a single 

106 coupling, using prederivatized Rink amide resin. Resin and all Fmoc amino acids were pur- 

107 chased from C S Bio Company Organic solvents and other reagents used for the sjmthesis 

108 and purification were high- performance liquid chromatography (HPLC) grade or better 

109 and purchased from Fisher Scientific (Ottawa, ON) and Sigma-Aldrich Canada (St. Louis, 

110 MO). Deprotection and cleavage of the peptides from the resin were conducted with a 

1 1 1 trifluoroacetic acid (TFA)/water (95:5 by voliraie) cleavage cocktail followed by cold precip- 

112 itation with tert-butyl ether. The crude products were piurified by preparative reverse-phase 

113 HPLC in a Vydac C-8 column by use of a water/acetonitrile linear gradient with 0.1% TFA 

114 as the ion pairing agent. The molecular weights of the peptides were confirmed by matrix- 

115 assisted laser desorption ionization time of flight (MALDl-TOF) mass spectrometry. The 

116 purified peptides were lyophilized and stored at 4 °C. 

117 Ala-Gly peptide with Mw=146.14 g/mole, Tyr-Leu peptide with Mw=294.35 g/mole, and 

118 SDS (99% purity) with Mw=288.38 g/mole were purchased from Sigma-Aldrich Canada 

119 (St. Louis, MO) and were used as received without further purification. Deuterium oxide 

120 D2O with 99.9% isotopic purity was purchased from Cambridge Isotope Laboratories (St. 

121 Leonard, Quebec). 

Table 1 Sample nomenclature. All samples were made with D2O as a solvent, and unless stated have 0.1 M 
sodiimn, oxalate buffer in them. Final concentrations [SDS] were achieved by mixing different stock solutions. The 
molar ratio R =[peptide]/[SDS]=30 was kept constant for GAD-2 solutions. 



Abbreviation 


Final [SDS] 


SDS-buf 

GAD-2-SDS 

Ala-Gly-SDS 

Tyr-Leu-SDS 


2-187 mM 

1- 80 mM 

2- 60 mM 

2-60 mM 
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1 22 G AD-2-SDS, Ala-Gly-SDS, Tyr-Leu-SDS, and SDS samples were prepared with compo- 

123 sitions according to table [T] The molar ratio (R) of SDS concentration to peptide concentra- 

124 tion in GAD-2-SDS samples was held constant (R = [SDS]/[GAD-2]=30). The concentration 

125 of dipeptides (Ala-Gly and Tyr-Leu) in Ala-Gly-SDS and Tyr-Leu-SDS systems was 2 mM. 

126 The pH value for all samples was adjusted to be 4 by the addition of sodium deuteroxide 

127 or deuterium chloride. All samples were made with D2O as solvent and, unless otherwise 

128 stated, have 0.1 M sodium oxalate buffer (Na2C204) in them. Sodium oxalate buffer was 



129 used in previous NMR studies to adjust the pH of SDS micelle-peptide solutions ( Orfi et al 



1998 Deaton et al 2001| . It is effective as a buffer for pH below 5, where the histidine-rich 



131 GAD-2 peptide is expected to have a net positive charge. Moreover, the chemical structure 

132 of sodium oxalate does not include protons in it. As a result, the one dimensional proton 

133 NMR spectra do not include buffer peaks that might overlap with SDS and peptides peaks. 

134 The self-diffusion measurements were carried out in a diffusion probe (DiffSO) and with 

135 maximum field gradient (1800 G/cm) at a resonance frequency of 600 MHz on a Bruker 

136 Avance II spectrometer. Diffusion was measured with a pulsed-field gradient stunulated- 

137 echo sequence ( Price 1997| with (almost square) trapezoidal gradient pulses. The diffusion 



138 coefficient of a molecule in aqueous solution is obtained from the attenuation of the signal 

139 according to the equation ( [Price 1997) 



l-|lil = -Dk (1) 



^S(O), 

140 where S(k) is the "intensity'of the signal (the integration of the relevant peak region) 

141 in the presence of field gradient pulse, S(0) is the intensity of the signal in the absence 

142 of field gradient pulse, k = (y6g)^(A-6/3) is a generalized gradient strength parameter, 

143 y = = 2.6571 X 10^ T"^s"^ is the gyromagnetic ratio of the nucleus, 6 = 2 ms is the 

144 duration of the field gradient pulse, A = 100 ms is the time period between the two field 

145 gradient pulses, and g is the amplitude of the field gradient pulse. Ti relaxation mea- 

146 surements were also performed using the DiffSO diffusion probe and the relaxation rates 



147 (Ri = 1/Ti) are reported in section 2.1 
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2 Results and Discussion 

Complementary NMR-based techniques were utilized in order to identify components for 
different samples based on their one-dimensional NMR spectra and to extract parameters 
such as self-diffusion coefficients and longitudinal relaxation rates. The one-dimensional 
(ID) proton NMR spectra at a resonance frequency of 600 MHz on a Bruker Avance II 
spectrometer and at sample temperature 298 K are shown in figure [T| In all cases the trace 
signal of HDO in D2O is the most dominant peak (at ~ 4.7 ppm); however the HDO, peptide 
and SDS peaks are all spectrally separable. 



12000 
10000 

8000 
6000 
4000 



10 8 6 4 2 [ppm] 



HDO 
GAD-2 




(a) 



(b) 



Fig. 1 ID iR NMR spectrum for (a) a peptide-free SDS sample with [SDS] = 6 mM (b) a GAD-2-SDS sample with 
[SDS]=60 mM and [GAD-2]=2 mM. Sample temperature is 298 K. 



In this work, we carried out experiments with peptide at varying SDS concentrations 
in the presence of sodium oxalate buffer. We also performed experiments on pure SDS 
solutions as well as buffered SDS solutions for comparison. Figure |2] shows the signal at- 
tenuation and the self -diffusion coefficients for SDS and peptides in a buffered peptide-free 
SDS sample and GAD-2-SDS sample. The signal attenuation in all samples was observed 
to be monoexponential. This suggests that the exchange of SDS molecules between the SDS 
in micelles and in free solution must be very rapid in the NMR time scale. The values of 
the observed diffusion coefficients were calculated from the monoexponential decays using 
equationjlj For peptide-free SDS solutions prepared with sodium oxalate buffer (figure[2a|. 
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Fig. 2 The attenuation of the signal S(k)/S(0) on a log scale versus k = (y6g)2(A- 6/3) for (a) a peptide-free SDS 
sample with [SDS] = 40 mM and 0.1 M sodium oxalate buffer (b) a GAD-2-SDS sample with [SDS]=80 mM, 
[GAD-2]=2.67 mM, and 0.1 M sodium oxalate buffer. 6 = 2 ms and A = 100 ms. 



165 the signal attenuation of SDS was obtained by integrating the area under the spectral region 

166 between to 4 ppm. For the GAD-2-SDS system, the spectral ranges from to 4 ppm and 7 

167 to 9 ppm were used to obtain SDS and GAD-2 signal attenuation, respectively. In each case 

168 the SDS and peptide spectral regions were chosen to ensure a clear spectral separation. 



169 2.1 Relaxometry 



1.4 
1.2 

£ 1.0 
of 

0.8 
0.6 

(a) 



4 i g 



SDS-buf 
fit 



Fit: yO+A*Exp(-[SDS]/C*) 
C*= (1.9±0.1)ml\l 
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Fig. 3 Proton longitudinal relaxation rates Ri for SDS in (a) peptide-free SDS solutions versus SDS concentration 
[SDS]. (b) GAD-2-SDS solutions with R=[SDS]/[GAD-2]=30. Similar measurements were carried out for the two 
dipeptides (not shown). The error bars represent the uncertainty in the numerical single exponential fit for the 
raw data of the relaxation time measurements. 
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170 Figure |3] shows the variation in the proton longitudinal relaxation rate Ri for the SDS 

171 peaks in peptide-free SDS solutions and GAD-2-SDS solution with increasing SDS con- 

172 centration. In all cases, the SDS concentration dependence can be fit to an exponential 

173 yO + A exp(-[SDS]/C*); where C* is identified as a characteristic SDS concentration that 

174 indicates a change in the local environment for the surfactant molecules in GAD-2-SDS 

175 and peptide-free systems. 

176 The exponential fit yields a C* - (1.9 ± 0.1) mM for peptide-free SDS samples (figure[3a|. 

177 For GAD-2-SDS samples, the spectral region from to 4 ppm is used to calculate Ri for the 

178 SDS peaks at different SDS concentration. This yields a C* = (1.3 + 0.1) mM (figure ??). 



179 2.2 Diffusometry 

180 2.2.2 Surfactant Solutions and Analysis Methods 



181 Figure 4a shows the self-diffusion coefficient of SDS in 3 peptide-free SDS systems: one with 

182 sodium oxalate buffer with pH=4 (red open circles), and two without sodium oxalate buffer. 

183 Of the unbuffered solutions one was with pH unadjusted but measured to be between 3 

184 and 3.5 (blue open squares), and one with the pH=4 (black filled squares). Below [SDS] 

185 = 60 mM, the SDS diffusion coefficient D^^ for unbuffered solutions at different pH are 

186 indistinguishable from each other, while values in the buffered solution are much lower. 

187 The pulsed-field-gradient signal attenuation is monoexponential for all samples (fig- 

188 ure|2j. This implies that the exchange of SDS molecules between the SDS in micelles and 



189 in free solution is rapid in the NMR time scale. Previous studies (Soderman and Stilbs 



1994[ | showed that a surfactant molecule visits more than one environment over very short 



191 timescales, and interpreted the observed diffusion coefficients using a two-site exchange 

192 model. In all the systems considered here, the SDS surfactant can either be a free monomer 

193 in solution or associated with a surfactant-rich aggregate. The observed self-diffusion coef- 

194 ficient of SDS is thus a linear combination of the self -diffusion coefficient D^^^ of the free 

tree 

195 molecules in bulk and that of the bound molecules in the micelle D^'^'^ „ in peptide-free 

micelle ^ ^ 

196 solutions or a peptide-SDS complex D^°^ 
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DSbs-Df°f, [SDS]<Co, 

^Ts = (D^°e - Di°y (fs) + DfgS^, [SDS] > Co (2) 

197 where fg = [SDS]free/[SDS] is the fraction of free SDS molecules, D^g^^. is either the micellar 

198 diffusion coefficient in peptide-free samples, or the diffusion coefficient of the micelle- 

199 peptide complex, and Co refers to the critical (micellar or aggregation) concentration (CMC 

200 or CAC), and [SDS] is the total SDS concentration. A key assumption of the model is that 

201 there are only two distinct species, the free and the aggregate states; however, as will be 
seen later, we are able to check for self -consistency of the model. 




0.1 1 10 100 1000 

[SDS] (mM) [SDS] (mM) 

(a) (b) 

Fig. 4 Self-diffusion coefficient in peptide-free SDS solutions, (a) D versus SDS concentration [SDS] for solutions 
with sodium oxalate buffer (pH=4) (red open circles), and unbuffered, with pH=3-3.5 (blue open squares), and 
with pH=4 (black filled squares), (b) Fraction (fg) of free SDS with and without sodium oxalate buffer. 



203 For simple spherical micelle systems, buffered and unbuffered peptide-free SDS solu- 

204 tions, [SDSJfj-ee = CMC for [SDS] > CMC. Therefore, equation|2]has 3 parameters, Co = CMC, 

205 D?°^ and D^'^^ „ . Fitting the buffered peptide-free SDS solution to the two-species model 

tree micelie o r r r 

206 in equation ^ yields the parameters D^ds = (4.90 ± 0.07) x 10"i° m^/s, D^°s = (6.3 ± 0.4) 

207 xlO"i^ m^/s, and CMC = (0.91 ± 0.02) mM, while for the unbuffered peptide-free SDS solu- 

208 tionD^DS = (4.71 + o.08) x 10-^° m^/s, D^^^^^ = (6.1 ± 0.9) x 10"" m^/s, and CMC = (5.3 ± 0.2) mM. 



Diffusion NMR Study of Complex Formation in Membrane- Associated Peptides 



11 



209 The main physical insight hidden in these curves is the onset of crowding. While the 

210 unbuffered and buffered solutions have very different djmamics at low [SDS], they converge 

211 to the same value above 60 mM signifying a concentration where hydrod5mamic interactions 

21 2 between complexes slow down the d5mamics and nullify the differences between aggregate 

213 ^size. Such an effect of hydrodynamic interactions has indeed been previously noted ( |Ando 

214 [and Skolnickl|20T0). 



215 2.2.2 Peptide: GAD-2 

When the size of a hydrophobic peptide is large enough that surfactant motion is rapid 
on the timescale of peptide motion, the peptide is expected to be associated with several 
surfactant molecules and there should never be free peptide, i.e. the peptide binding fraction 
is close to 1 . For example, in the G AD-2-SDS system, since the concentration of SDS is 30 
times higher than GAD-2 concentration (R = [SDS] /[GAD -2]- 30), we know that there is 
unlikely to be free peptide. In this case, D^'^g^ = D'^^P'"^'^. Using this additional information 
allows us to use the two-site model even if the D^^ versus 1/[SDS] relationship is not 
linear. The only proviso is that the overall particulate volume fraction must always be small 
enough that hydrodynamic effects are negligible. For the peptide-free SDS system, we have 
seen that this is generally true for concentrations below 60 mM, or volume fractions below 
0.04. For GAD-2-SDS system, the size of an GAD-2-SDS aggregate is expected to change 
with SDS concentration. Therefore, the concentration [SDS]free of free SDS monomers is 
expected to change in the SDS concentration regime above CAC. We may simply rewrite 
and rearrange equation|2]for [SDS] > Co but with DSDS^^D^^P'i'i^ 

rSDSlf DSps_DPep«de 

free 



216 Figure 5a shows the self -diffusion coefficient of GAD-2 and SDS in the GAD-2-SDS 

217 system. We measured the self-diffusion of GAD-2 in the SDS concentration range that is 

218 higher than 13.3 mM. Due to experimental limitations (small value of signal-to-noise ratio). 
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GAD-2-SDS 
O SDS 
■ GAD-2 

— two species model 



; D^J° =(3.6 ± 0.2) X 10"" m';s I 
. Df,°,^ =(5.0±0.2)x10"'°m';s 
CAC=(0.73 ± 0.03) mM 
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Fig. 5 (a) Self-diffusion coefficient of GAD-2 and SDS in a GAD-2-SDS system witfi R=[SDS]/[GAD-2]=30 versus 
SDS concentration [SDS] (b) Fraction (fs) of free SDS and concentration ([SDS]free) of free SDS versus SDS 
concentration [SDS] 



219 we were not able to extract tlie self -diffusion coefficient of GAD-2 in tlie SDS concentration 

220 range below 13.3 mM, but we were able to measure the surfactant diffusion. 

221 The SDS self diffusion coefficient is fit well to the two species model for [SDS]<25 



222 mM (figure 5a solid line), and it deviates from the fit for higher SDS concentration (fig- 
dotted line). The resulting fit parameters are D|°^ = (5.0 ± 0.2) X 10''^'^ m^/s, D|DS 



223 ure 



5a 



224 = (3.6 ± 0.2) X 10"^^ m^/s, and CAC= (0.73± 0.03) mM. 



225 PFG-NMR can be used to get spectrally-resolved diffusion coefficients. Using both SDS 

226 and peptide diffusion coefficients as a function of [SDS], we extract the fraction (fg) of free 

227 surfactant in the monomer state in the aqueous solution as well as the concentration of free 



228 surfactant [SDS]fi.ee- This is shown in figure 5b With increasing surfactant concentration, 

229 fs decreases wliile [SDS]f].ee rises from 0.7 mM (close to the CAC) to ~ 1 mM (close to 

230 the CMC). This is consistent with the picture ( [Barhoum and Yethiraj 2010 Jones 2002 1 

231 that the concentration of free surfactant above the CAC/CMC is equal to the value of the 

232 CAC/CMC . In the peptide-SDS system, and similar to the behavior in the nonionic polymer- 



anionic surfactant (poly(ethylene)oxide-SDS) system (Barhoum and Yethiraj 2010 1, the free 



234 concentration rises further until it reaches the CMC value in the buffered solution. 
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Next, we estimate the hydrodynamic radius Rh of GAD-2-SDS complexes (figure 6a I 
using the Sutherland - Stokes - Einstein equation (Jones 2002) 

235 where Kb is the Boltzmann constant, T is the absolute temperature, and r| is the solvent 

236 viscosity (rio^o-l-l mPa.s). The hydrod5mamic radius Rh is obtained from the peptide 

237 diffusion (D = dP^P'^'^^ open squares in figure 



6a I as well as from the fitted value of D^'-'^ 

Aggr 



238 obtained from the concentration dependence of the surfactant diffusion (dashed red line 

239 in figure [6a). For [SDS] < 25 mM the hydrodynamic radii obtained via peptide diffusion 

240 and surfactant diffusion are roughly the same, with a value of approximately 5.5 nm. 

241 Interestingly, Rh (obtained from peptide diffusion D^'^P*^'*'^) increases as a function of SDS 

242 concentration to about 10 nm at 60 mM, less than a factor of two increase. 



243 Plotted in figure 6b is the variation in the ratio of SDS molecules to peptide molecules in 

244 a complex r - ([SDS] - [SDS]free)/([SDS]/R) - (1 -fg) R exhibits a very slight increase, from 

245 ~ 28 to 29, and approaches R = 30 asymptotically. We need to imderstand how the aggregate 

246 size changes in order to accommodate the two-fold increase in the hydrodynamic radius 

247 Rh; we will return to this point later. 



248 2.2.3 Comparison with smaller dipeptides 

249 In order to study the effect of peptide size on the dynamics of peptide-SDS complexes, 

250 diffusometry was carried out to quantify complex formation of SDS with the dipeptides 

251 Ala-Gly and Tyr-Leu. 

252 A plot of the SDS self-diffusion coefficient for all systems in the current study in one 

253 graph (figure |7j shows clearly that SDS diffusion looks similar for the systems with small 

254 di-peptides (Ala-Gly and Tyr-Leu) and the peptide-free SDS system with sodium oxalate 

255 buffer This suggests that the fraction of free SDS in the Tyr-Leu-SDS and Ala-Gly-SDS 

256 systems is similar to the fraction of free SDS in the peptide-free SDS system with buffer 

257 (figure On the other hand, SDS diffusion looks very different for the system with 
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Fig. 6 (a) The apparent hydrodynamic radius (Rj^ ), extracted from the peptide diffusion coefficient, of 
GAD-2-SDS complexes versus SDS concentration [SDS]. The horizontal dashed line is the value of apparent 
hydrodynamic radius (R^'''' =5.5 ± 0.3 nm) obtained via the SDS aggregate diffusion coefficient 
(D|°|^ = (3.6 ± 0.2) X 10~i^ m^/s ) at low SDS concentration [SDS]< 25 mM, vi^here the two species model is valid. 
This value is in agreement with the peptide diffusion coefficient. At high SDS concentration [SDS] > 60 mM. The 
diffusion coefficient measured gives no information about the true hydrodynamic radius. The intermediate SDS 
concentration regime, denoted by the gray area, is the regime in which either complex size is indeed increasing 
with concentration or hydrodynamic interactions between complexes is slowing down the motions, (b) The ratio 
of SDS molecules to peptide molecules in a complex (r) versus SDS concentration ([SDS]) for GAD-2-SDS 
samples. 
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Fig. 7 Self -diffusion coefficient of SDS in all the systems studied: peptide-free SDS system with sodium oxalate 
buffer, Tyr-Leu-SDS, Ala-Gly-SDS, and GAD-2-SDS with R=[SDS]/[GAD-2]=30 samples. 



258 long peptide (GAD-2-SDS system), suggesting that the GAD-2-SDS complexes are very 

259 different from the Ala-Gly-SDS and Tyr-Leu-SDS complexes, which are essentially indis- 

260 tinguishable from micellar aggregates with no peptide. This means that the peptide-micelle 

261 binding characteristics of the Tyr-Leu and Ala-Gly dipeptides are different from the much 

262 longer GAD-2 peptide. Also, this indicates that GAD-2 significantly disrupts the micellar 

263 aggregate. This conclusion likely extends to other long and hydrophobic peptides. 
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264 3 Conclusion 



265 NMR-based techniques have been utilized in this work to study the nature of peptide- 

266 micelle complexes in a buffered 19-residue antimicrobial peptide (the GAD-2-SDS system). 

267 First, we examined the impact of the buffer ( figure [4a|. Varying the pH over a small range 

268 in the absence of a buffer shows no effect on the micellar structure, while the CMC is lower 

269 in the presence of the buffer. 

270 At high surfactant concentrations ([SDS] > 60 mM), the observed diffusion coefficients 

271 of surfactant SDS molecules merge for buffered and unbuffered peptide-free SDS systems. 

272 In addition, the linear two species model (equation |2j, is valid below [SDS] =60 mM but 

273 not valid above. This is similar to the findings in previous work for a system of anionic 

274 surfactant (SDS)-nonionic polymer polyethylene oxide (PEO) ( [Barhoum and Yethiraj 2010 1 

275 where this concentration was identified as the onset of macromolecular crowding: this 

276 refers to the excluded volume effect of one macromolecule with respect to another ( [Zhou 



et al 2008). Our primary finding is that [SDS] =60 mM signals the concentration beyond 



278 which one cannot, even in principle, extract hydrodynamic radii or aggregate ratios. 

279 At low surfactant concentrations ([SDS] < 25 mM), the observed diffusion coefficient of 

280 SDS (figure [5a] is well described by the two-species model in equation [2| Moreover, in this 

281 range, the surfactant aggregate diffusion coefficient and the peptide diffusion coefficient 

282 coincide. This is a self-consistency check that gives confidence in the linear two species 

283 model and the results obtained. 

284 At intermediate SDS concentrations, the apparent hydrodjmamic size increases from 5.5 



285 nm at 25 mM to 10 nm at 60 mM (figure 6a I. This increase in the apparent hydrodynamic 

286 size might either reflect a true increase in aggregate size, or it might indicate the existence 

287 of hydrod5mamic interactions between complexes. Given that the ratio of SDS to GAD- 

288 2 molecules in a complex is not changing by much, i.e. r ^ R (figure |6b|, an increase in 

289 the mean aggregate size might arise from an increase in the average number of peptides 

290 in one complex from 1 (at 25 mM) to approximately 2 (at 60 mM). A third possibility is 

291 that such an increase in hydrod5mamic radius arises from a change in shape (for example 
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292 from spherical to oblate or prolate) i jBloomfield 2000| . However, in order to account for 

293 a factor two increase, one would need to have a rather spectacular shape change with 

294 a formation of very anisotropic complexes with an approximately 20:1 axial ratio. These 

295 three possibilities - an increase in number of peptides in a complex, long-range interactions 

296 between complexes, or a dramatic change in complex shape - are depicted in figure |8] 

297 As noted by (Zhou et al 2008 Schreiber et al 2009}, a deeper understanding of role of 

298 electrostatic and hydrod5mamic interactions is needed in the study of macromolecular 

299 crowding, and this needs to be studied further. 



o o 



oo c 

O CD 



or 




[SDS] < 25 mM 25 mM < [SDS] < 60 mM 

Fig. 8 A schematic diagram showing each peptide-surfactant complex as a single isolated complex (left, isolated 
circles) at low SDS concentrations. Results at intermediate SDS concentrations are consistent with either two 
peptides in each complex schematically represented by two circles (top right), highly anisotropic complexes 
(right middle), or long-range hydrodjmamic interactions represented by arrows between complexes (right 
bottom). 

300 There is a distinct difference in the behavior of large peptides (M^'^^'^'^'^ > M^''^^'^'^"') and 

301 small dipeptides (M^'^P'''*'' ~ M^^^^*^"') ■ The small dipeptides ( Ala-Gly and Tyr-Leu) hardly 

302 affect the SDS diffusion coefficient (figure [7||. This indicates that the dipeptides behave just 

303 as the surfactant does: i.e. rapidly exchanging between aggregate and free state. For large 

304 peptides such as GAD-2, on the other hand, rapid exchange between free and aggregate 
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305 state is practically impossible. This is because the surfactant molecules form micellar-like 



306 aggregates along the peptide chain, consistent with a bead-on-a-charn picture ( Chari et al 



2004 Roscigno et al |2003| for large-molecule aggregates. We therefore expect the approach 



308 outlined in this work to be valid generally for large hydrophobic peptides. 

309 As a conclusion, some recommendations are suggested in order to extract parameters 

310 such as the hydrodynamic radius and peptide binding fractions. All our results consistently 

31 1 show that measurements should be made well below the concentration that signals the onset 

312 of crowding (about 60 mM for SDS). Furthermore, even at intermediate concentrations, the 

313 true hydrod5mamic size of the peptide-SDS complex is not necessarily constant, and the 

314 concentration dependence of the hydrodynamic radius can still not be ignored. The only 

315 unambiguous concentration-independent statements can be made at low concentrations: 

316 in this system, this is below [SDS]= 25 mM. 
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